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Spinless fermions on highly frustrated lattices are characterized by a lowest 
single-particle band which is completely flat. Concrete realizations are provided 
by the sawtooth chain and the kagome lattice. For these models a real-space 
picture is given in terms of localized states. Furthermore, we find a finite zero- 
temperature entropy for a suitable choice of the chemical potential. The entropy 
is computed numerically at finite temperature and one observes a strong cooling 
effect during adiabatic changes of the chemical potential. We argue that the 
localized states, the associated zero-temperature entropy and thus also the large 
temperature variations carry over to the repulsive Hubbard model. The relation 
to flat-band ferromagnetism is also discussed briefly. 
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1. Introduction 

Frustrated quantum magnets exhibit a rich variety of semi-classically ordered and 
disordered ground states. This has been analyzed in some detail for two-dimensional 
models (see [1,2] for recent reviews). Interesting behaviour is also observed in a 
magnetic field, including plateaux in the zero-temperature magnetization curve and 
field-induced quantum phase transitions (compare [2-4] and references therein). A 
particularly intriguing phenomenon are non-interacting localized magnons which 
have recently been discovered at the saturation field of the spin-s XX Z model on 
highly frustrated lattices [2-6]. This gives rise to an enhanced magnetocaloric ef- 
fect. Indeed, cooling rates around the saturation field in geometrically frustrated 
classical spin systems were predicted [7] to be up to several orders of magnitude 
bigger than in non-frustrated spin models. Enhanced cooling rates have also been 
verified experimentally, namely for GdsGasOi^, a hyper-kagome lattice (see [8] and 
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references therein), and the pyrochlore magnet Gd2Ti 2 07 [9]. Since these results 
suggest applications for efficient low-temperature magnetic refrigeration, the mag- 
netocaloric properties of the corresponding quantum systems are currently under 
intense investigation [10-15]. 

With the aid of the Jordan- Wigner transformation (see e.g. Chapter 1 of [16]), 
s — 1/2 XXZ chains can be mapped to spinless fermions. These are generally in- 
teracting, but the interaction is irrelevant at low densities for sufficiently dispersive 
bands such that the low-temperature magnetocaloric properties at the transition to 
saturation in some quantum spin chains can be understood in terms of free spin- 
less fermions [10]. Although there is no such direct connection for flat bands or in 
higher dimensions, we may still expect that free spinless fermions capture relevant 
qualitative features. This motivates us to analyze low-energy properties of free spin- 
less fermions on the sawtooth chain (section |3J) and the kagome lattice (section EJ), 
after a brief summary of some basic equations in section EJ In section El we discuss 
implications for the Hubbard model on the same lattices and mention a relation to 
flat-band ferromagnetism [17-21]. 

2. Model and basic equations 

A central quantity of this paper is the entropy S. It is given in terms of the 
partition function Z by 1 



Furthermore, we will analyze free spinless fermions, whose Hamiltonian reads 



The first sum runs over the nearest-neighbor pairs of a lattice with N sites, 
cj and di are fermion creation and annihilation operators at site i with anticommu- 
tation relations {c|,Cj} = S it j, and the number operator at site i. The 

hopping parameter £y > corresponds to the exchange constant of an antiferro- 
magnetic spin-s XXZ model, and the chemical potential \i to a magnetic field h 
along the z-axis (the sign of /i has been chosen positive in order to allow for a direct 
comparison with the XXZ model in a magnetic field). We denote the single-fermion 
energies of (0) by and note that this one-particle problem is equivalent to the one- 
magnon problem relative to the ferromagnetically polarized state of the spin-s XX 
model on the same lattice. Furthermore, the Hilbert space of spinless fermions on 
an iV-site lattice is isomorphic to that of the s = 1/2 XXZ model on the same 
lattice. Nevertheless, the many-particle states are not equivalent. Apart from the 
different statistics, interactions are absent in the model of spinless fermions (J2J). The 
latter property makes it straightforward to derive the thermodynamics, which is a 
challenging problem for the XXZ model. 

lr The Boltzmann constant will be set to unity throughout this paper fc^ = 1. 
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Figure 1. The sawtooth chain. The hopping parameter along the baseline is t, 
while along the diagonal directions it is t' . For t' = \[2t the one-particle states 
can be localized, as indicated by the bold line and the coefficients of the wave 
function ((HJ) written next to the three associated sites. 



For a free Fermi-gas the grand partition function 2 is given by Z = Yl k (l + e £fe / T J 
Substitution into (0) yields the entropy for a system of free fermions 

Sm =2>(l + e-/>iE^££. (3) 

k 1 k 1 ^ ^ 



3. Sawtooth chain 

The first model which we will discuss is the sawtooth chain shown in Fig. ^ The 
one-particle problem for (J2J) on the sawtooth chain leads to a 2 x 2 matrix whose 
eigenvalues yield two bands of one-particle energies 

e±(Jfe) =t cos(k)±^t 2 cos 2 0) + 2t' 2 cos(k) + 2t' 2 + p . (4) 

When one chooses hopping parameters along the diagonals t' = v2 1 in terms of the 
hopping parameters along the baseline t, the one-particle energies become 

e-(k) = -2t + n, e+(k) = It (1 + cos(A;)) (5) 

Note that the lowest band £-(k) is completely flat. The result (jSJ is equivalent to 
the one-magnon energies of an XX spin model on the sawtooth lattice [3]. Very 
similar results are also obtained for the Heisenberg model on this lattice (compare 
eq. (33) of [10] for an explicit expression of the one-magnon energies for s = 1/2). 

Now we will concentrate on the case t' = \[2t where the lowest band £-{k) is 
completely flat. As in the XX Z model [3,4], it is possible to localize these one- 
particle states on the three sites indicated by the bold line (a ,, valley") in Fig. [T] 
To be more precise, let us number the sites along the baseline by 2i and the ones 
at the top of the triangles by 2 i + 1 . If we further denote the fermionic vacuum by 
|0), the localized wave functions are given by (up to normalization) 

\i 2i ) = 4. |o> , lI = 4_! + 4 +1 - V2 4 . (6) 

At // = 2t we have £-(k) = 0. This implies a degeneracy of the ground states with 
< n < 1/2 due to the 2^ distinct ways to fill the zero-energy band. In this 

2 The chemical potential /i is included in the definition of the one-particle energy Sk for state k. 
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Figure 2. Lines of constant entropy for free spinless fermions on the sawtooth 
chain with t' = \[2t. The value of the entropy per site S(T)/N is indicated next 
to each line. 



This corresponds to the zero-temperature entropy found in the Heisenberg model 
at the saturation field. However, the value in the latter case is smaller, namely 
S(0)/N = 0.2406059. . . [10-12]. That a larger value is obtained for free fermions 
can be understood as follows. Many-particle states can be constructed from lo- 
calized one-particle states in different valleys. For spinless fermions it is possible 
to put one-particle states in neighbouring valleys. The corresponding two-particle 
state |0) is an exact eigenstate with energy 2e_. However, in the XX Z 

model the energy of two localized magnons occupying neighbouring valleys is higher 
since they share a site [10,12]. This condition that two localized magnons may not 
occupy two neighboring valleys lowers the ground state degeneracy and thus the 
zero-temperature entropy of the XX Z model. 

The entropy at finite temperature and with arbitrary [i can be evaluated from (jSJ). 
One can either go to the thermodynamic limit and replace the sums by integrals. 
Alternatively, one can choose a large system size (N = 2048, for example) and 
evaluate the sums numerically. The latter has been performed to obtain Fig.|2]which 
is indistinguishable from the thermodynamic limit iV — > 00. The right half of this 
figure can be compared to Fig. 6(a) and Fig. 7 of [10] which show numerical results 
for the entropy of the s = 1/2 Heisenberg model in the /i-T-plane. The following 
different low-temperature regions can be distinguished in Fig. |2J (i) for \i > 2 1 both 



manner we have identified a zero-temperature entropy at [i — 2t, t' = y2t 



SO n2 

= = 0.34657... . 

N 2 



(7) 
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Figure 3. The kagome lattice. The bold hexagon indicates a localized one-particle 
state, the signs in the corresponding wave function are written next to the six sites 
of this hexagon. Circles indicate the closest non-overlapping packing of localized 
one-particle states. 



bands are empty and all excitations are gapped, (ii) for \i = 2t the ground state is 
degenerate with the T = entropy given in (JJJ), (iii) for 2t > fi > the chemical 
potential is in the band gap between the lower band £-(k) and the upper band 
£+(k) such that excitations again have a gap, (iv) for > fi > —At the chemical 
potential is in the upper band e + (k) and one has gapless low-energy excitations, (v) 
for n < — A t both bands are filled and all excitations are again gapped. 

Now let us start at a small temperature T and change the chemical potential 
adiabatically. Such an adiabatic process keeps the entropy constant, therefore the 
temperature T has to change along curves as in Fig. |21 If one approaches region 
(iv) either from region (iii) or region (v) one sees that temperature drops as fi —>■ 
or fi — > —At. An even larger temperature drop is observed when one approaches 
point (ii) either from region (i) or region (iii). This enhanced cooling is clearly due 
to the zero-temperature entropy J7J) at [l = 2t. In fact, if the entropy at the starting 
point is below this value, one theoretically approaches T — > as fi — > 2t. The 
qualitative behavior is similar to that of the XX Z model on the sawtooth chain [10] 
(an interacting model) despite quantitative differences. 



4. Kagome lattice 

In this section we discuss a two-dimensional lattice, the kagome lattice shown in 
Fig. El Since a unit cell of this lattice contains three sites, the one-particle problem 
for the Hamiltonian (0) now leads to a 3 x 3 matrix whose eigenvalues yield three 
bands of one-particle energies 

e {k x ,k y ) = /I -2t, 



£±(k x ,k y ) = fj, + t±t 



l + 4co S (|)(cos(4^) + cos(|)). ,8) 
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Figure 4. The three bands £i(k) of single-particle energies for spinless fermions 
on the kagome lattice along the path in the Brillouin zone shown in the inset. 
Note that £o(k) is completely independent of k. 



Fig. |3] shows these three one-particle bands along a path in the Brillouin zone. Note 
that the lowest band e (k) turns out to be completely flat. These results can be 
related to known ones. Firstly, using a particle-hole transformation [17], the present 
problem maps to a tight-binding problem of electrons moving on the kagome lattice 
(compare Fig. EJe.g. with Fig. 1(b) of [22]). The one-magnon problem for the XX Z 
model on the kagome lattice is another equivalent problem (see [3] and in particular 
Fig. 2.31 of [2], compare also eq. (1) of [12]). 

As in the XX Z model [2,3], it is again possible to localize these one-particle 
states in real space, namely on a hexagon, as indicated by the bold line in Fig. El 
If we number the sites around a hexagon h with j = 1, . . . , 6, these localized wave 
functions can be written as (compare also Fig. |3]for the signs (— l) 3 ') 



\h) = 4 |o> 



(9) 



The energy of the lowest band becomes s (k) — at \i = 2t. Consideration of all 
possible ways to fill this zero-energy band yields 2 N ' ?J zero-energy excitations. For 
\i — 2t another one-particle energy vanishes, namely £_(0, 0) = (see Fig. 0}. How- 
ever, this occurs only at a single point k = 0. Therefore, the additional contribution 
can be neglected in the thermodynamic limit N — > oo. Hence, the zero-temperature 
entropy per site is at // = 2 1 

S(0) In 2 

^-^ = — = 0.231049.... (10) 

The XX Z model on the kagome lattice also gives rise to a finite zero-temperature 
entropy at the saturation field [2,11,12]. A mapping of the configurations Q to the 
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Figure 5. Lines of constant entropy for free spinless fermions on the kagome 
lattice. The value of the entropy per site S(T)/N is indicated next to each line. 

exactly solved problem of hard hexagons [23,24] yields a lower bound S(0)/N > 
0.11108 in this case [11,12]. On the one hand, one expects to find a lower entropy in 
the XX Z model than for spinless fermions, because creation of localized magnons on 
two neighbouring hexagons costs no energy for free spinless fermions while localized 
magnons can be created for the XX Z model only on spatially separated hexagons. 
On the other hand, the situation remains more complicated than in one dimension 
because of additional topological effects. In fact, a class of many-magnon states has 
been found in the XX Z model which are independent of hard-hexagon configura- 
tions [12,14]. Therefore we believe that the precise value of the zero-temperature 
entropy at the saturation field of the s — 1/2 Heisenberg model on the kagome 
lattice deserves further attention. 

Also for the kagome lattice it is straightforward to obtain the entropy at finite 
temperature and with arbitrary [i from © and the single-particle energies (jHJ). We 
have again evaluated the sums numerically on a large but finite lattice (N = 10800). 
The corresponding values of the entropy per site are shown in Fig. and can be 
considered as representative for the thermodynamic limit. The following four low- 
temperature regions can be distinguished in Fig. El (i) fc> r t 1 > 2 t all bands are empty 
and all excitations are gapped, (ii) for \i = 2t the ground state is degenerate with 
the T = entropy given in (fTU|) . (iii) for 2t > // > — At the chemical potential is in 
one of the bands e±{k) and one has gapless low-energy excitations, (iv) for ji < —At 
all bands are filled and all excitations are again gapped. Additional structures in 
region (iii) can be attributed to features in the T = single-particle density of states 
(compare Fig. Note that in contrast to the sawtooth chain, the kagome lattice 
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does not give rise to a gap for \i slightly below 2 t. 

If one now changes \i adiabatically, temperature drops considerably as fi — > 2 t. 
One can even reach T = during an adiabatic process if the starting entropy is less 
than the T = entropy (jl(Jj) at fi = It. Note the asymmetry of the curves in Fig. |5] 
around fi = 2t: The dependence on /i is generally weaker in the gapless regime (iii) 
than in the gapped regime (i). One also observes a cooling effect as fi — » — At from 
the gapped regime (iv), although unlike in Fig. El there is no pronounced minimum 
close to /i = — At. 

5. Hubbard model and flat-band ferromagnetism 

Now we will discuss the influence of interactions using the Hubbard model as 
example. The fermions are assigned an additional spin index a =f, |, and there is 
an on-site Coulomb repulsion U > for fermions (electrons) with different spin: 

N N 
o" i=l o" i=l 

The Coulomb repulsion is not effective for spin-polarized configurations (all a equal). 
These configuration can be identified with spinless fermions. Conversely, each eigen- 
state of the spinless fermion Hamiltonian (J2J) yields a spin-polarized eigenstate of 
the Hubbard model (jlljl . In particular the one-particle problem is identical in both 
models. The single-particle bands of the Hubbard model are simply twofold degener- 
ate. Using results from the previous sections, we can now establish a T = entropy 
for the Hubbard model (jllj) for U > and otherwise the same conditions which 
lead to T = entropy for spinless fermions (0). A relation to flat-band ferromag- 
netism [17-21] is also evident. 

Let us assume that the lowest single-particle band is completely flat, as we found 
for the sawtooth lattice (Fig. [TJ) at t' = y/2t and the kagome lattice (Fig. BJ). The 
corresponding states can be localized in real space as before; the creation operator 
of a local state lX a just acquires an additional spin index. States localized at ji, j 2 , 
. . . with spins a±, cr 2 , . . . can be created by repeated application of distinct creation 
operators 

4,,w- - () )- (12) 

These states are clearly eigenstates for any ji ^ j 2 if they are spin-polarized aj = a. 
As long as the local states are created in spatially non-overlapping regions, the 
states ()12j) with <j\ ^ cr 2 are also eigenstates with the same energy as the spin- 
polarized states. However, once the localization regions overlap, there are generally 
two electrons with distinct spin at a site. The Coulomb repulsion in (fTTj) raises the 
energy of such states by an amount oc U. Note that certain combinations where local 
states with different spins overlap must have the same energy as the spin-polarized 
states because of the 5'?7(2)-symmetry of the Hamiltonian (fTTj) . However, if Coulomb 
repulsion is active, it will raise the energy of many-particle states in the lowest band. 
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Now choose n such that the lowest band has zero energy e(k) = (// = 2 1 for the 
sawtooth chain and the kagome lattice). The preceding discussion then yields a lower 
bound of the T = entropy of the Hubbard model (fTTj) in terms of spinless fermions 
©: SHubbard(O) > Ssf(O). The precise value of the entropy of the Hubbard model 
(jllj) is difficult to determine because of the additional zero-energy states which are 
not spin-polarized. In any case, also in the Hubbard model adiabatic changes of \i 
will give rise to large temperature variations. 

Spatial overlaps of the individual local states cannot be avoided at large filling 
fractions of the flat band. Spin-polarized states still have the lowest possible energy in 
this case while a non-polarized state is generally pushed higher in energy by Coulomb 
repulsion. This implies a phenomenon known as flat-band ferromagnetism [17-21]: In 
this region of densities spin-polarized, i.e. ferromagnetic configurations are preferred 
over non-magnetic configurations. To be more precise, let us consider the kagome 
lattice. Using our conventions (fTTj) . the lowest band can be filled in a spin-polarized 
manner for < n < 1/3. On the other hand, local wave functions have to overlap for 
n > 1/9 (see Fig. |HJ). Slightly above the threshold n = 1/9, ferromagnetism can still 
be avoided by having overlap only for localized states with equal spins. According 
to a detailed mathematical analysis [17], the ground state of the Hubbard model 
(fTTj) is ferromagnetic for 1/6 < n < 1/3. 

6. Summary and Outlook 

In sections |3] and 0] we have examined spinless fermions on the sawtooth chain 
and the kagome lattice. We found that the lowest single-particle band is completely 
flat for a particular choice of hopping parameters in the case of the sawtooth chain 
and without any fine-tuning for the kagome lattice. These states can be localized in 
real space. Since the particles have zero energy at // = 2 1, one finds a finite T = 
entropy at this value of the chemical potential. This implies large adiabatic tem- 
perature variations when [i is changed adiabatically at finite temperature. At least 
on the single-particle level, the construction is completely analogous to the XXZ 
model on the same lattice. It can therefore be generalized to many other lattices 
for which localized magnons have been constructed [2-4,6], including in particular 
the pyrochlore lattice in three dimensions. We do not wish to speculate about pos- 
sible experimental applications. Rather we regard these non-interacting models as 
an illustration of important qualitative features of spin models. First experimental 
investigations [8,9] of the magnetocaloric effect in highly frustrated magnets shows 
that they are indeed promising candidates for efficient low-temperature magnetic 
cooling. 

Tuning of the chemical potential /i through a non-degenerate band-edge, i.e. 
an edge without flat directions, yields a second-order quantum phase transition. 
Spinless fermions thus provide a simple realization of the quantum critical scenario 
discussed in detail in [25,26]. 

In section we have briefly commented on the Hubbard model where an on-site 
Coulomb repulsion U gives rise to interactions. However, this repulsion is not active 
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for spin-polarized configurations. Therefore, the existence of a finite T = entropy 
for a certain value of \i and consequently the large temperature variations when \i is 
tuned through this point carry over from spinless fermions to the Hubbard model. 
Nevertheless, both the T = entropy and the cooling rate in the Hubbard model 
remain to be analyzed in more detail. In this context it is also intriguing to note 
that a finite T = entropy arises in the Hubbard model under the same conditions 
which lead to flat-band ferromagnetism [17-21]. 
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